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Comment on "Critical State Theory for Non- 
parallel Flux Line Lattices in Type-II Supercon- 
ductors" 

Recently, Badia and Lopez (BL) proposed a new model 
relevant to the electrodynamics of hard superconductors 
in multicomponent situation [fij. BL have claimed that 
using their least action principle with the given constrains 
on the current density one can solve any critical state 
problem. They argue that the BL model is equivalent to 
the set of Maxwell equations with the vertical current- 
voltage characteristics (CVC). However, for the isotropic 
case (the only case analyzed in the BL formalism 
does not lead to new results in comparison to the Bean 
critical state model generalized for the multicomponent 
case. As regards the anisotropic situation, we show here 
that the BL model is not equivalent to the set of Maxwell 
equations with the "vertical CVC". In other words, we 
call in question the validity of the BL model in a general 
case. 

BL proposed a new formalism for solving the criti- 
cal state problems and applied it for the analysis of the 
isotropic situation. For this particular case, the BL for- 
malism is surely adequate and is completely equivalent to 
the set of Maxwell equations with the CVC in the form, 
J = J C E/E, where J c is the critical current density, E 
is the electric field. This material equation corresponds 
to the generalized Bean critical state model (GBCSM). 
The GBCSM is commonly used in the electrodynam- 
ics of hard superconductors. Specifically, the analysis of 
this model allowed us to discover a new interesting phe- 
nomenon of the collapse of the transport current S and 
of the magnetization pi . Contrary to [pj , the GBCSM is 
useful for clarification the nature of the collapse. Indeed, 
let us consider a dissipation-free dc critical current Jdc 
shielding external dc magnetic field H and providing dc 
magnetization M of the superconductor. After switch- 
ing on ac magnetic field h _L H, ac current J ac is induced 
(J ac -L Jdc)- According to GBCSM, the total current 
should flow along E. This means that the dissipation- 
free current Jdc is suppressed and the collapse of M ap- 
pears to be. It is natural that owing to equivalence of 
the BL and GBCSM models, simulations concerning the 
collapse in jlj] should coincide with the results obtained 
previously (see e. g. B). 

Although Letter jjj deals with an analysis of the 
isotropic case, BL declare that "several extensions may 
be implemented if dictated by the physics of the prob- 
lem. These include the effect of equilibrium magneti- 
zation by means of an appropriate B(H) relation, the 
selection of the model J C {H) and the use of anisotropic 
control spaces" . It seems to us that the validity of the 
BL formalism in the anisotropic case is questionable. Al- 
though the BL model operates correctly with the region 
of the dissipation-free current flow, we argue that the 
set of Maxwell equations with CVC is not equivalent to 
the BL model for the dissipative region. Indeed, CVC of 


samples in the critical state contains not only informa- 
tion about a region where the current flow is dissipation- 
free, but provides the one-to-one correspondence between 
the directions of the critical current and the electric field 
within the dissipative region. It should be noted that not 
only a single CVC but a class of CVC can correspond to 
the same restriction on the current density. For example, 
the class of CVC with an arbitrary parameter a, 

Jx = J c cos(9 + a), J y = J c sin(# + a), (1) 

corresponds to the isotropic constraint |J| < J c but de- 
scribes the case where the direction of J is inclined with 
respect to E. Here, x and y are arbitrary coordinate axes 
in the sample plane, 8 is the angle between the vector E 
and x axis. Obviously, the set of Maxwell equations to- 
gether with these CVC would give different profiles for 
the field penetration, whereas the BL formalism with the 
isotropic constraint gives a unique solution of the prob- 
lem. This means, that the variational procedure chooses 
only one of CVC from a class, ignoring the real CVC. 
This choice is defined by the requirement of the mini- 
mum of the BL functional. Only for the isotropic case 
such a requirement predetermines a = in Eq. ([!]) (i. e. 
J||-E) owing to the symmetry of the problem. In other 
cases, one cannot be sure that the suggested least action 
principle chooses the CVC correctly. Of course, if this 
principle is proved, it would be a very important tool 
to determine the CVC via the known restriction region. 
Unfortunately, there is no reason to expect this because 
the minimum principle for the entropy production rate is 
not proved for nonlinear unsteady systems. 

Thus, the BL formalism being valid for the isotropic 
case can hardly be correct for more complex situations, 
in particular, for anisotropic cases. 
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